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ABSTRACT 

A detailed analysis of XMM-Newton observations of C1J0046. 3+8530 {z = 
0.624) is presented. The cluster has a moderate temperature {kT = 4.1 ± 
0.3 keV) and appears to be relaxed. Emission is detected at > ia significance 
to a radius of 88% of i?200 (the radius enclosing an overdensity of 200 times 
the critical density at z = 0.624) in a surface-brightness profile. A temperature 
profile (extending to 0.7i?2oo)j and hardness-ratio map show no significant de- 
partures from the global temperature, within the statistical limits of the data. 
The cluster's bolometric X-ray luminosity is Lx = (4.3 ± 0.3) x 10"*^ erg s~^, 
and the total mass derived within -R200: assuming isothermality and hydro- 
static equilibrium, is A/200 = 3.0^°^ ^ IO^^Mq. The gas entropy measured 
at 0.1i?2O0 is compared with a sample of local systems, and found to be con- 
sistent with self-similar evolution with redshift. The metallicity, gas density 
profile slope, and gas mass fraction are all consistent with those found in low-z 
clusters. 

Key words: cosmology: observations - galaxies: clusters: general - galaxies: 
high-redshift - galaxies: clusters: individual: (C1J0046. 3-1-8530) - intergalactic 
medium - X-rays: galaxies 
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1 INTRODUCTION 

The X-ray properties of galaxy clusters, particularly at 
high redshift, are very sensitive probes of cosmology 
(e.g. Henry 1997; Jones ct al. 1998; Arnaud ct al. 2002; 
lAllen et alJl200a) . Cluster mass is the property most 
directly related to the predictions of cosmological mod- 
els, but also one of the most difficult to measure. X- 
ray observations directly measure the properties of the 
gaseous intra-cluster medium (ICM), which can be used 
to derive the total gravitating mass of a cluster, pro- 
vided that the gas is in hydrostatic equilibrium. This 
derivation also requires knowledge of the temperature 
structure of the gas. Thus to measure accurate cluster 
masses with X-ray observations, detailed measurements 
of the dynamical state of the ICM and its temperature 
structure are needed. Such observations of local clus- 
ters are becoming common, but detailed studies, and 
particularly measurements of ICM tem perature struc- 
ture at high z (z > 0.6) , are very rare 1 Jeltema et alJ 
I2OOII: lArna-ud et ai]l200l iMauehan et alJl2004^ . 
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C1J0046. 3-1-8530 was discovered as an extended X- 
ray source in the Wide Angle RO SAT Poin ted Survey 
fWARPSi lScharf et al.lll997: Pcrlman et al.;2002') . and 
confirmed as a galaxy cluster with redshift z = 0.624. 
The identification was based on its extended X-ray na- 
ture, an overdensity of faint galaxies, and concordant 
optical spectroscopic redshifts of two galaxies, includ- 
ing the the galaxy coincident with the X-ray centroid 
JHorner et al.ll2004r . The redshift of the cluster is con- 
firmed below via an X-ray redshift measurement. Prior 
to the observation detailed here, relatively little was 
known about the cluster's properties. C1J0046. 3+8530 
falls 11' off-axis in two consecutive XMM-Newton obser- 
vations of the open star cluster NGC 188, taken on the 
29-30 October 2000 (observation identifiers 0100640101 
and 0100640201) using the medium filter. The analysis 
of these data are described below, and the properties of 
the system are derived and discussed. 

A ACDM cosmology of Hq = 70 km s"^ Mpc"^ , 
and SIm = 0.3 (JIa = 0.7) is adopted throughout, and all 
errors are quoted at the 68% level. At the cluster's red- 
shift, 1" corresponds to 6.8 kpc in this cosmology. The 
virial radius is defined here as -R200, the radius within 
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which the mean density is 200 times the critical density 
at the redshift of observation. 



2 OPTICAL OBSERVATIONS 

The first optical follow-up observation of 
C1J0046. 3+8530 was obtained in August 1997 in 
the form of a 3 x 4 minute exposure in the R band with 
the University of Hawaii's 2.2m telescope. Mediocre 
atmospheric conditions combined with a high airmass 
of 2.2 resulted in poor seeing of 1.2 arcsec as measured 
off the combined final image. Deeper imaging was 
performed subsequently in August 2001 in the R and 
I bands (40 and 20 min respectively) with the 4.2m 
William-Herschel Telescope (WHT). Again, the seeing 
was 1.2 arcsec. All the imaging data were reduced 
in the s tandard manner . Observations of 6 standard 
stars of iLandold (|l992) were used to calibrate the 
WHT photometry. Galaxy number counts made in 
an area excluding the cluster w ere c onsistent with 
those of iKQmmel fc Waeneil (1200 J) and iMetcalfe et alJ 
(jlSSS)! confirming the accuracy of the photometry, and 
indicated that the limiting magnitude was R = 22. 

Spectroscopic observations of galaxies in the field 
of C1J0046. 3+8530 were performed on the Keck 10m 
telescope using t he LRIS low-dispersion spectrograph 
ijOke et alJ 1199,4 1 in August 1997. The instrumental 
setup combined a 2048x2048 Tektronix CCD with 
24^m pixels (LRIS-R CCD2), the 300/5000 grating, the 
GG495 order blocking filter, and a longslit of 1.5" width 
to obtain a dispersion of 2.5 A/pixel over a spectral 
range extending from 5150 to 8950A. Two galaxies were 
observed for a total of 900 seconds. The airmass of the 
targets was 2.4, the seeing 1.0 arcsec. Standard image 
processing techniques were used to reduce the obtained 
data, and the measured redshifts are given in Table Q 



3 DATA PREPARATION 

The two XMM-Newton observations of C1J0046.3+8530 
were filtered, cleaned and prepared in the standard 
way (see [Maughan ct al. (2004) for details), to give two 
events lists for each detector. Both observations were 
affected by brief periods of high background, and their 
removal left 30 ks (observation 1) and 50 ks (obser- 
vation 2) of good time. After cleaning, the mean PN 
count rates in the 10 — 15 keV energy band (excluding 
bright sources) we re consistent with that found in the 
iLumb et al.l ll2002ll blank-sky background files, indicat- 
ing that the high-background periods were successfully 
removed. 

The position of C1J0046.3+8530 near the edge of 
the field of view raised important calibration issues. The 
90% encircled energy radius of the Point Spread Func- 
tion (PSF) at 1.5 keV increases from 47" on-axis to 52" 
at the cluster centroid for the MOS detectors, and from 
51" to 65" for the PN detector. These values were de- 
rived with the Calview tool in the XMM-Newton Sci- 
ence Analysis Software (SAS). Following recent advice 
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Figure 1. The radial dependence of the PN effective area at 
1.5 keV. The dashed line marks the cluster position. 



from the XMM-Newton help-desk^ these encircled ener- 
gies were computed using the 'EXTENDED' accuracy 
setting in Calview. This is the most accurate method for 
estimating encircled energies, and is valid at all off-axis 
angles. These encircled energy radii are therefore more 
up-to-date than the values given in the XMM-Newton 
user's handbook v2.1. In addition, the shape of the PSF 
changes from being radially symmetric at the optical 
axis, to a more arc-like shape, elongated orthogonally 
to the radial direction. In our modelling of the cluster's 
surface-brightness distribution, where the shape of the 
PSF is important, model PSFs are computed with the 
'MEDIUM' accuracy setting in Calview, which provides 
the only two-dimensional models of the PSF. 

The effective area also decreases strongly with off- 
axis angle. Fig.^shows the dependence of effective area 
on off-axis angle; at 700", the off-axis angle of the clus- 
ter position, the effective area has fallen to less than 
50% of the on-axis value. These instrumental issues are 
accounted for as much as possible in the following anal- 
ysis, and their possible effects are discussed. 

In addition, the energy-dependent vignetting func- 
tion of the telescopes has recently been verified by ob- 
servations of compact superno va remnants at a ra nge of 
off-axis angles and azimuths JLumb et al.M2004f) . This 
analysis revealed a misalignment of « 1' in the opti- 
cal axis of the telescopes containing the PN and MOS2 
cameras. For the cluster location in this observation, the 
possible discrepancy in the assumed vignetting would 
mean that the PN fiux has been overestimated by ft; 6% 
and the MOS2 flux underestimated by ~ 8%. A sec- 
ondary effect of the energy dependence would be that 
the corresponding temperatures would also be slightly 
over- and under-estimated respectively. These effects are 
not accounted for in the following analysis, as the uncer- 
tainties introduced are small compared to the statistical 
uncertainties. 



^ See| http://xnim.viIspa.esa.es/xmnihelp/Calibration?id=9356| 
for a discussion of PSF issues. 
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Figure 2. An exposure-corrected image of the XMM-Newton 
field of view. Data from both observations are combined, and 
binned with image pixels of 4.4". The inset is a blow-up of 
the dashed region, and clearly shows CIJ0046. 3-1-8530 and 
XMMU J004751. 7-1-852722 as extended sources. 




85:30:00 



0:48:00 47:30 0:47:00 46:30 0:40:00 45:30 0:45:00 

Right ascension 

Figure 3. Contours of X-ray emission detected by XMM- 
Newton (0.3 — 5 keV) overlaid on combined University of 
Hawaii 2.2 m and WHT R-band images of C1J0046. 3-1-8530. 
The contours are logarithmically spaced above the lowest 
contour at 0.32 counts pixel^^, which is 1.5 times the back- 
ground level. The embedded panel is an enlargement of a 
central 1' X 1' region, showing the positions of the galaxies 
with measured redshifts 



4 X-RAY IMAGING ANALYSIS 

Images and exposure maps were produced for the PN 
and MOS cameras for both observations in the en- 
ergy range 0.3 — 5 keV to maximise signal-to-noise. The 
exposure maps for each instrument were "spectrally 
weighted" by creating several exposure maps produced 
in narrow energy bands, in which the energy dependence 
of the vignetting function was small. These maps were 
weighted by the photon flux predicted in each energy 
band by the best-fitting spectral model to the global 
cluster emission (see Q and summed. The images and 
weighted exposure maps for each camera were then com- 
bined to give a mosaiced image and map. The expo- 
sure map was normalised so that its value at the clus- 
ter centroid was 1, and the image was divided by the 
result. This removes the effects of vignetting and chip 
gaps, while maintaining as far as possible the Poisso- 
nian statistics of the data. An exposure-corrected image 
of the field of view, combining the data from both ob- 
servations is shown in Fig|5| The diffuse emission from 
ClJ0046.3-f8530 is visible near the Northern edge of the 
field of view, and the increasing PSF distortion with off- 
ajds angle is apparent. 

The exposure-corrected image was adaptively 
smoothed so that all remaining features were signifi- 
cant at > 99% confidence, and contours of the smoothed 
emission were overlaid on an R-band optical image, as 
shown in Fig. |3 Ignoring the distortion of the contours 
caused by the two point sources to the South, the con- 
tours appear fairly circular, indicating that the cluster's 
atmosphere is reasonably relaxed. 

The two galaxies with measured redshifts are la- 
belled in Fig. 121 and their positions and redshifts are 
given in Tabled The table also gives the position of the 
X-ray centroid and the measured X-ray redshift of the 



cluster (see i|SJ. All of the following assumes z — 0.624, 
the redshift of the central galaxy, as the cluster redshift. 

4.1 An extended source to the South East 

A region of extended emission (designated XMMU 
J004751. 7-1-852722) was detected 4' South East of 
C1J0046.3-I-8530 (see Fig.EJ. Contours of this emission 
(smoothed to 99% significance) are overlaid on an op- 
tical image in Fig. 0] The centroid of this extended 
emission coincides with the 2MASS galaxy 2MASX 
J00475079+8527213 (K = 13.5 ± 0.2 mag). Within 
a radius of 32" from the source centroid, a total of 
~ 500 net photons were detected, using a local back- 
ground annulus, and excluding point sources and emis- 
sion from ClJ0046.3-f8530. The properties of XMMU 
J004751. 7+852722 are discussed in 37|and ^ 



4.2 Modelhng the surface-brightness 
distribution of C1J0046. 3-1-8530 

The surface-brightness distribution of C1J0046. 3-1-8530 
was modelled in Sherpa with a two-dimensional (2D) 
method. A 0.3 — 5 keV image from each detector (with- 
out exposure correction) was fit with a surface bright- 
ness model which was multiplied by an appropriate ex- 
posure map, and convolved with a PSF. Due to the low 
nu mbers of cou nts in a typical pixel, the Cash statis- 
tic iCashlll97Sr) was used when fitting the models. The 
cluster emission model used was a 2D paramctcrisation 
of the /3- model (Cavalicrc & Fusco-Fcmiano 1976). Two 
additive background model components were used. The 
first was a simple constant-background model which, 
after multiplication by the exposure map should ac- 
curately represent a flat, but vignetted X-ray back- 
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Object 



» [2000.0] 



<5 [2000.0] 



Galaxy A 

Galaxy B 

C1J0046. 3+8530 



00^46"" 17. 73= 
00'^46""25.92= 
00'i46™17.83= 



+85°31'13.6" 
+85°30'34.1" 
+85°31'09.5" 



0.6243 ± 0.0004 
0.630 ±0.002 



0.615 



+0.008 



Table 1. Positions and redshifts of the two galaxies marked in Fig. Inland of the cluster itself. The position given for the cluster 
is the X-ray centroid, and the redshift is the X-ray redshift (derived in JSJ- 
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Figure 4. Contours of a smoothed X-ray image of XMMU 
J004751. 7+852722 from both observations and all EPIC 
detectors overlaid on a WHT R-band image. The con- 
tours are logarithmically spaced above the lowest contour at 
0.25 counts pixel"^, which is 1.5 times the local background 
level. 



ground. The second background component consisted 
of a 0.3 — 5 keV image produced from closed-filter data, 
to account for the particle background. This particle 
image was divided by the exposure map, so that after 
multiplication by the exposure map in the fitting pro- 
cess, the background model cons isted of a vignetted and 
a non- vignetted component (see iMaughan et al.l (|2004l l 
for a discussion of the background components and these 
models) . 

Due to the large wings of the PSF at the detector 
position of the cluster, it was not possible to fully ex- 
clude point source emission without excluding almost all 
of the cluster emission. Instead, the point sources within 
the fitting region were modelled as delta functions, and 
their position and amplitude were fit to the data. This 
procedure is not ideal, because the fitting software only 
allows the model being fit to be convolved with one PSF 
per dataset, while the shape of the PSF varies over the 
detector region spanned by the cluster emission (as can 
be seen in Fig. EJ. Calview was used to generate PSF 
images for the detector position of the cluster centroid 
at an energy of 1.5 keV for each detector (the PSF is 
not strongly energy dependent in the 0.3 — 5 keV band). 
The 'MEDIUM' accuracy level was used, which provides 
the only available 2D model of the XMM-Newton PSF, 
including its arc-like distortion at the cluster off-axis 
angle. 

These models were fit simultaneously to the PN and 
MOS images from each observation. The results of the 
fits were not consistent between the observations. The 



best-fitting /3-profile core radius {tc) and slope (/3) are 
-+^**" and /3 = 0.49l^:g| in observation 1, and 
and P = 0.60to;o3 in observation 2. This 
discrepancy was investigated in detail, and is likely to be 
caused by a combination of the large PSF, uncertainties 
in its modelling, and a variation in the flux of bright- 
est point source between the observations. The count 
rate from the bright point source to the South East was 
25% higher in the first observation. Given the uncer- 
tainties in modelling the PSF, it is likely that there is 
more point source contamination remaining in observa- 
tion 1, increasing the number of counts at large radii, 
which would cause the /3-profile slope to fit to a lower 
value. In addition, the core radii found in both observa- 
tions are small compared to the PSF, which means they 
may not be reliably measured. This would also affect 
the values of /3 as there is degeneracy between the two 
parameters. These fits were repeated with the core radii 
fixed at 22" , which corres ponds to 150 kpc. a val ue typ- 
ical of local clusters (e.g. Ijones fc Formanlll999f) . With 
this constraint, the best-fitting values of all other pa- 
rameters are consistent between the observations. 



With this in mind, we chose to use only observa- 
tion 2 in our imaging analysis because it was the longer 
observation, with better signal-to-noise, and the point 
source contamination is lower. The best-fitting model 
has re = 20.2l:ts" and /3 = O.eoto'og, with a low el- 
lipticity of 0.07, and we proceed with these values for 
all further analysis. With r^ fixed at this value, fitting 
to the observation 1 data gives /3 = 0.57 ± 0.05, consis- 
tent with the above result. Fig. |S] shows a radial profile 
of the best-fitting 2D model, including instrumental ef- 
fects, plotted on a profile of the data. The point sources 
were not excluded from the profile of the data in this 
comparison, because they are included in the 2D model. 
The dotted line is a profile of the best-fitting background 
model components. These profiles show that the model 
is a reasonably good description of the data, although 
it is clear again that the point sources at ~ 50" and 
~ 100" are not perfectly modelled. The dip in the pro- 
files at X 130" is due to the edge of the PN CCD. This 
is apparent in both the data and the models because 
the models have been multiplied by exposure maps to 
replicate the instrumental effects in the data. 

In order to measure the extent of the detection, 
a one-dimensional radial profile of the X-ray emission 
was extracted from the combined (both observations, 
all cameras) exposure corrected image. Point sources 
were excluded out to a radius of 20", and the profile 
was adaptively binned so that all radial bins were sig- 
nificant at the 3(7 level above the background. It is not 
appropriate to fit a one-dimensional /3-model to this pro- 
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Figure 5. A radial profile of the best-fitting 2D model to the 
C1J0046. 3+8530 observation 2 data, including instrumental 
effects, is shown with the solid line. The points are a radial 
profile of the data, and the dotted line is a profile of the 
best- fitting background model components. 
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Figure 6. MeKaL temperature with corresponding hard- 
ness ratio values for different hard and soft energy bands, 
derived from simulated EPIC spectra. The Galactic ab- 
sorption, redshift and instrumental responses appropriate to 
C1J0046.3+8530 were used. 



file, because this would not account for the large PSF. 
Emission is detected at the 3o" level out to a radius of 
142" (965 kpc), which corresponds to 88% of -R200, as 
discussed below (©. While some point source contami- 
nation will remain after their exclusion (due to the large 
PSF wings), the diffuse emission is clearly extended to 
~ 140" at azimuths with no point sources. 



4.3 Hardness-Ratio Map 

The temperature structure of the cluster was probed 
with a hardness-ratio map. The hardness ratio, HR, was 
defined as 



HR = 



H - AH, 



bg 



S - ASbn 



(1) 



where H and S are the counts in the source region in 
the hard and soft band respectively, and the bg subscript 
indicates the counts found in a background region. A is 
the ratio of the area of the source region to the back- 
ground region. The error on the hardness ratio is then 
given by 



a{HRY 



H + A'^Hbg {H^AHtgfiS + A^Si 

{s-ASt^r^^ 



(S - AStgV 



^ 



Befor e defining the hard and soft bands, absorbed 
MeKaL JKaastra &: Mewd llQQST spectra were simu- 
lated for the three EPIC detectors (with the appropri- 
ate redshift. Galactic absorption and instrumental re- 
sponse) to estimate the sensitivity of HR to gas tem- 
perature variations, with different hard and soft energy 
bands. A soft band of 0.3 — 1 keV, and a hard band of 
1 — 5 keV were chosen because HR values in the range 
0.5 — 2 (within which the statistics are not dominated by 
one band with low counts) corresponded to a tempera- 
ture range of « 1 — 10 keV, a reasonable range about the 
global temperature. Fig. |S| illustrates this, showing the 
relationship between hardness ratio and MeKaL tem- 
perature for different hard and soft bands. 

An adaptive binning template was computed from 
a broad-band image, to define image bins with > 
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Figure 7. Adaptively binned HR significance map of 
ClJ0046.3-f 8530 data. The overlaid contours are the same as 
in Fig. 13 and the scale-bar shows the HR significance colour 
scaling. 



100 counts in order to minimise statistical uncertain- 
ties while keeping reasonable spatial resolution. A few 
bins were allowed to fall below this threshold to im- 
prove the resolution, but all bins with < 50 counts 
were excluded from further analysis. The background- 
subtracted, exposure-corrected hard- and soft-band im- 
ages were binned according to this template, and hard- 
ness ratios with errors were computed for each bin. Fig. 
|7|shows a significance map where the value in each bin is 
given by the difference in HR between the data at that 
point, and the HR corresponding to the global temper- 
ature, divided by the error on the HR in that bin. 

All of the bins (excluding those associated with the 
point sources) are consistent within the statistical un- 
certainties with the emission being isothermal; 17 of 25 
(68%) are within la of the global HR and all bins are 
within 2a. However, as the larger image bins in the 
HR significance map are similar in size to the 75% en- 
circled energy radius of the PSF, this method would 
only be sensitive to reasonably strong temperature vari- 
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ations. Of the two point sources, that to the South West 
(a[2000.0] = 00''45'"53.8P, (5[2000.0] = +85°30'30.6") 
has harder emission (at the 2.7a level), which may 
be consistent with it being an absorbed AGN in the 
foreground galaxy that appears to be the counterpart 
of the X-ray point source in Fig. The South East 
point source (a[2000.0] = 00''46™57.62", (5[2000.0] = 
+85°30'36.6") has significantly (~ lOcr) softer emission, 
and appears to be associated with a foreground star. 
The spectra of these sources and the galaxy cluster are 
discussed in JS] 



5 X-RAY SPECTRAL ANALYSIS 

Spectra were extracted for all EPIC detectors for both 
observations. For comparison, background spectra were 
extracted both locally, and from blank-sky background 
files. The results of the spectral analysis were found to 
be independent of the background type used, so local 
background spectra were used in all further analysis, as 
the systematic uncertainties involved were considered 
to be smaller. The source spectra were extracted within 
a circle of radius 88" centred on the cluster centroid, 
corresponding approximately to the 3ct detection radius 
in each separate observation (as opposed to the Scr de- 
tection radius of 142" when the observations were com- 
bined). The strong radial dependence of effective area 
at the cluster position (Fig. was taken into consider- 
ation when selecting the regions from which to extract 
the local background spectra. Exposure maps of each 
detector were used to select spatial regions where the 
mean effective exposure (equivalent to effective area) 
was close to that in the source spectrum extraction re- 
gion. Ancillary response files (ARFs), weighted by the 
spatial distribution of counts in the source region, were 
also generated for the cluster's off-axis position. 

The two point sources to the South fell within the 
source region, and their emission was excluded to a ra- 
dius of 20" (~ 75% encircled energy radius of the PSF 
of all XMM-Newton telescopes at 0.3 - 5 keV). Due to 
the large PSF wings, and its elongation at this off-axis 
angle, some point source contamination will remain, and 
the effect of this is examined below. 

The spectra from the three EPIC detectors in 
the 0.3 — 5 keV band were each grouped to give > 
20 counts per energy bin, and were fit simultane- 
ously with an absorbed MeKaL model. This was done 
for the two observations independently, and the best- 
fitting model parameters were consistent. All spectra 
for both observations were then fit simultaneously with 
an absorbed MeKaL model. The best-fitting temper- 
ature and metallicity are kT — 4.1 ± 0.3 keV and 
Z = 0.39lgi3 respectively, when the absorbing col- 
umn was frozen at the Gal actic value of 7.6 x 10^" cm~^ 
JPickev fc LockmanllToOOl) . with x^/i' = 372/360. This 
model is shown with the data in Fig. |H| When the ab- 
sorbing column is allowed to vary, the best-fitting value 
is (6.8 ± 1.2) X 10^" cm~^, consistent with the Galactic 
value. In addition, an X-ray redshift was measured, and 
the best-fitting value is 2 = 0.615tQQo6' which supports 




channel energy (keV) 

Figure 8. All EPIC spectra for observations 1 and 2 
of C1J0046. 3-1-8530, and the best-fitting absorbed MeKaL 
model. Absorption is frozen at the Galactic value. For pre- 
sentation purposes, the data from the same instruments in 
the two observations have been combined for this plot, and 
adjacent bins have been combined to give a minimum signal- 
to-noise ratio of 3. The lower panel shows the ratio of data 
to model values. 



the optical redshift of 0.624 ascribed to the cluster based 
on the two measured galaxy redshifts. 

The spectrum of the soft, bright point source to 
the South East was examined, using a local annulus as 
a background region, to ensure that the cluster emis- 
sion is also included in the background spectrum. The 
spectrum is well fit by an unabsorbed MeKaL model 
ix^/f = 9.95/10) of kT = 0.53lg;i;^ keV and Z = 
0.06 ± 0.03Zq. The redshift of this model was frozen 
at zero. A power-law model is rejected at the 95% level 
by the data, supporting the conclusion that the emission 
is stellar. The unabsorbed flux from this point source is 
(8.4±1.9) X lO"^'' erg s"^ cm"^ in the 0.5-2 keV band. 

To test whether the emission from this star affected 
the cluster spectrum, a cluster spectrum was extracted 
including this source, and its emission was included in 
the model. The model used was an absorbed MeKaL 
model at z = 0.624 plus an unabsorbed MeKaL model 
at z — 0. The best-fitting temperatures of the cluster 
and star are kT = 4:.8tli keV and 0.38lj;;J^ keV re- 
spectively, with x^/'^ ~ 112.4/121. The good agreement 
between these fits, and those found when fitting the 
sources independently indicates that the cluster emis- 
sion is not seriously contaminated by the star, and 
most importantly, that the temperature measured for 
the cluster is not influenced by any remaining stellar 
emission. 

A cluster spectrum was also extracted with the 
emission from the harder South- West point source in- 
cluded. When this spectrum was flt with an absorbed 
MeKaL model, the best-fltting temperature was 4.9 ± 
0.5 keV, hotter than the best-fitting temperature when 
that point source was excluded. This is consistent with 
the harder spectrum of that source. When an additional 
absorbed power-law component with a spectral index of 
1.5 was included in the model, the best-fitting MeKaL 
temperature reduced to 3.8 ± 1.0 keV, which agrees well 
with the cluster temperature measured when the South- 
West point source was excluded. This suggests that the 
cluster spectrum was not significantly contaminated by 
South- West point source after it was excluded. 
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As a further test, the cluster emission model was 
fit to the data with the point sources excluded, in the 
energy range 2 — 8 keV. In this band the contribution 
from the soft stellar source should be negligible. The 
best-fitting temperature is kT = 4.11q g keV, providing 
further evidence that the measured temperature was not 
affected by point source contamination. 

The flux of C1J0046. 3+8530 measured by ROSAT 
in the 0.5 — 2 keV pa ssband is f .6 x fO"^"^ erg s~^ cm~^ 
iJHorner et alJ 120041) when extrapolated to an infinite 
radius. For comparison, the XMM-Newton flux in this 
band is (f .45 ± 0.07) x fO"" erg s"^ cm"^ (after ex- 
trapolation to infinite radius). This difference in flux is 
due to the two point sources, which were unresolved (or 
barely resolved in one case) in the ROSAT data, and 
included in the .Horner ct al.. (.2004') flux measurement. 



5.1 Temperature Profile 

With '-^ 5000 net counts (all EPIC detectors combined) 
within the detection radius, it was possible to measure a 
temperature profile of the gas in C1J0046. 3-1-8530. Spec- 
tra were extracted in 3 annular bins for all detectors, 
and both observations, and fit simultaneously with an 
absorbed MeKaL model, with the abundance frozen at 
O.3Z0, and the absorbing column frozen at the Galac- 
tic value. A locally extracted background spectrum was 
used. The best-fitting temperatures are shown in Fig. 
|5| and are consistent (with large errors) with the global 
temperature of 4.1±0.3 keV, though the outer bin could 
indicate cooler gas at larger radii. The spectral fit to 
the data from the outer bin was also performed in the 
2 — 8 keV band, to test the effect of any contaminating 
emission from the bright point source. The best-fitting 
temperature agrees well (with large errors) with that 
shown in Fig. |5| The data do not suggest any cool gas 
in the central regions. 

The effect of emission from gas in the outer annuli 
being projected onto the inner annuli was modelled with 
an 'onion skin' method. The temperature structure was 
modelled as a series of spherical shells (each of which 
was isothermal), and the spectra were fit from the out- 
ermost shell in. The spectrum of a shell was modelled 
with a single temperature MeKaL component, plus a 
MeKaL component for each external shell, whose tem- 
perature was fixed at the value measured in that shell, 
and whose normalisation was multiplied by a scale fac- 
tor. These scale factors accounted for the volume of each 
external shell along the line of sight to the shell being fit, 
and the variation in density across each external shell 
using the measured gas density profile. The deprojected 
temperatures are also plotted in Fig.|^ The form of this 
deprojected profile is not significantly different from the 
projected profile, and does not reveal any central cool 
gas, although the size of the errors is increased, as one 
would expect, as there were less photons available to 
constrain the temperature of the free component in the 
interior bins. 

It should be noted that the effect of the large PSF 
mixing the emission between annular bins was not taken 
into account here. The 90% encircled energy radius of 
the PSF at the position of the cluster at 1.5 keV is ~ 65" 



projected k I ^ 
deprojected kT t 



Radius (arcsec) 

Figure 9. Projected and deprojected temperature profiles of 
CIJ0046. 3+8530, based on spectra fit with abundance frozen 
at O.SZq, and Galactic absorption. The solid line is the best- 
fitting global (emission weighted) temperature, with la er- 
rors represented by the dashed lines. The temperature error 
bars have been off'set by 2" on the deprojected data points 
for clarity. 



for the PN detector. This means that it is possible there 
was a temperature gradient that was blurred out by the 
PSF. We derived the total gravitating mass below as- 
suming isothermality, which is consistent with the ob- 
served deprojected temperature profile, but we also dis- 
cuss the possible effects of a temperature gradient. 

This temperature profile enabled the measurement 
of the entropy profile of the cluster gas, where en- 
tropy (S) is defined here as S = T jnj keV cm^. 
The electron density profile was derived from the best- 
fitting surface-brightness profile, normalised to the cen- 
tral value derived from the meas ured MeKaL normali- 
sation (see lMaughan et alJ i2003r ). and assuming spher- 
ical symmetry. The entropy profile is shown in Fig. 1101 
The entropy at O.I-R200 (with 7?200 defined as described 
in SlSJl is 180 ± 24 keV cm^, with the error taken from 
the fractional temperature error on the nearest bin in 
the deprojected profile. This is lower than the value 
of ~ 300 keV cm^ fou nd in local systems o f the same 
global temperature bv lPonman et all ll2003ll . This dif- 
ference in entropy is compared with the predicted self- 
similar evolution in f|5| 



6 DETERMINATION OF GLOBAL 
PROPERTIES 

The properties of C1J0046. 3+8 530 were derived follow- 
ing the methods described in iMauehan et al.l (|200^. 
In summary, the total mass profile was derived assum- 
ing hydrostatic equilibrium, spherical symmetry, and 
isothermality, and used to define an overdensity pro- 
file with respect to the critical density at the redshift of 
observation. The virial radius, -R200, was then defined 
as the radius corresponding to an overdensity of 200. 
The central gas density was computed from the mea- 
sured MeKaL normalisation, and used to normalise the 
gas density profile. The errors on all the non-measured 
quantities were derived by 10, 000 Monte-Carlo ran- 
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Figure 10. Entropy of C1J0046.3+8530 is plotted against 
radius as a fraction of -R200 ■ The solid line shows the entropy 
derived assuming isothermality at the global temperature, 
and the data points correspond to the dcprojccted tempera- 
ture bins in Fig. 1^ 
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Figure 11. Colour-magnitude diagram for C1J0046. 3-1-8530 
using galaxies within 60" of the X-ray centroid. Galaxy B 
in ClJ0046.3-f8530 with a spectrally confirmed redshift is 
marked as a solid square. The line shows the predicted red 
sequence a,t z = 0.624 (see text). 



domisations based on the errors on the measured prop- 
erties. 

Our assumption of isothermality is supported by 
the measured temperature profile and hardness ratio 
map, while the circular appearance of the X-ray con- 
tours and low measured ellipticity support the assump- 
tions of hydrostatic equilibrium and spherical symme- 
try. 

The derived mass profile of C1J0046. 3-1-8530 en- 
abled the definition of 7?200 ~ l-lOlyoy Mpc, which 
means that the cluster emission is detected out to 
0.88i?2oo and the temperature profile is measured out to 
0.687?2oo- This greatly reduces the uncertainties in ex- 
trapolating the cluster properties to i?200- The gas mass 
fraction of the system at i?200 is found to be 16 ± 4%. 
The properties of ClJ0046.3-f8530 are summarised in 
Table H 



7 OPTICAL IMAGING ANALYSIS 

The R and I band WHT images were used to investi- 
gate the properties of the galaxies in ClJ0046.3-f 8530 
and XMMU J004751. 7-1-852722. Galaxies were selected 
by using SExtractor (Bcrtin & Arnouts 1996;), and 
their R — I colours and total magnitudes were mea- 
sured. Colours were measured using apertures cho- 
sen to maximise signal to noise, and total magni- 
tudes were obtained from the SExtractor MAG_AUTO 
measurements. Magnitudes and colours were corrected 
for Galactic extinc tion using the measurements of 
ISchleeel et al.l il998r . The galaxies within a radius 65" 
of the X-ray centroid of the XMMU J004751. 7-^852722 
and within 60" (0.4i?2oo) of the X-ray centroid of 
ClJ0046.3-f8530 were then selected. These are plot- 
ted on a colour-magnitude diagram in Figs. 1111 and 
m Galaxy B in C1J0046.3-I-8530, with a spectrally 
confirmed redshift, is marked separately. The bright- 
est cluster galaxy (galaxy A in Fig. |3J was par- 
tially blended with a bright star, leading to large 
uncertainties on its colour and magnitude, so is not 
shown here. There is a possible red sequence at i? — 



I « 1.3 in the C1J0046. 3-1-8530 galaxies, while a 
stronger sequence a.t R — I ~ 1 is clear in the galax- 
ies of XMMU J004751. 7-^852722. This indicates that 
XMMU J004751. 7-1-852722 is at a lower redshift than 
C1J0046. 3-1-8530, perhaps unsurprising given the coin- 
cidence of the source with the relatively bright 2MASS 
galaxy. 

To estimate the redshift of XMMU 
J00475 1 . 7+85272 2, the stellar po pulation synthe- 
sis models of Bruzual A. fc Charlod (l993) were used 
to predict the normalisation of the red sequence vs. 
redshift. A single stellar population was used with a 
burst of star formation at z/=3 and passive evolution 
thereafter, a Salpeter IMF and solar metallicity. These 
are typical v alues found for cluster elliptical evolution 
studies {e.g. lEUis fc .Ton es 2004). Similar results are 
obtained if z/=2 is assumed. The predicted colours 
were normalised to th ose of the Coma cluster at L* 
lEisenhardt et alJ (|2fl04^ . The resulting predictions have 
been plotted in Figs. ^^ and ^| for redshifts z=0.624 
and z=0.25, using appropriate valu es for the slope 
from the Coma cluster compilation of lEisenhardt et alJ 
(200_^. The prediction matches the observed colours 
of the red sequence galaxies in C1J0046. 3-1-8530, and a 
redshift of z=0.25 is a good match to the red sequence 
in XMMU J004751.7-H852722. 

The BCG in XMMU J004751. 7-^852722 has 
R=17.4 and KRil3.5, at the limit of the 2MASS sen- 
sitivity. The R-K colour is consistent with that of local 
elliptical galaxies when the k-corrections are taken into 
account. Assuming the redshift of z = 0.25 is correct, 
the absolute magnitude is Mr ~ — 23.4-|-51og(/i7o), and 
this is a luminous giant elliptical galaxy, comparable in 
luminosity to BCGs in some rich, high Lx systems. 



8 PROPERTIES OF THE NEARBY 
CLUSTER XMMU J004751. 7+852722 

A hardness ratio was computed for XMMU 
J004751. 7+852722, using the energy bands defined 
above, giving HR = 1.74 ± 0.24. This HR value was 
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Cluster 



Redshift T(keV) Li„i(erg s"!) rd kpc) 



B.20o{Mpc) Mgas(R20o)(MQ>) Mtot(R20a)(MQ,) 



C1J0046.3+8530 



0.624 



4.1 ±0.3 (4.3 ± 0.3) X lO* 



1371 



0.601 



1.10;° ;°7 (4.7 ± 0.6) X lO'-'' 3.0ij; ;? x 10" 



Table 2. Summary of the measured and inferred properties of cluster C1J0046. 3+8530 based on XMM-Newton observations, 
assuming a cosmology of Qjvf = 0-3 {^a = 0.7) and Hq = 70 km s~^ Mpc~^ . The luminosity and masses were scaled out to 
^200- 



—it 



Figure 12. Colour-magnitude diagram for XMMU 
J004751. 7+852722 using galaxies within 65" of the X-ray 
centroid. The line shows the predicted red sequence at 
z = 0.25 (see text). 



converted to a temperature using simulated spectra as 
described above, with Galactic absorption and a metal- 
licity of O.3Z0. As the redshift of the source was not 
certain, spectra were simulated at a range of redshifts, 
giving a range of possible temperatures. At z = the 
HR value corresponds to kT — 3.0lj'5 keV, increasing 
to kT = i.Stl'l keV at the estimated redshift z = 0.25. 
The measured count rate was then converted to a 
0.5 — 2 keV unabsorbed flux, using these derived 

2 
2 



temperatures, giving Fx ~ 5.4 x 10 ^^ erg s ^ cm" 
{kT = 3.0 keV) and Fx = 6.2 x 10"^^ erg s"^ cm" 
{kT = 4.8 keV). The corresponding bolometric X-ray 
luminosity at z=0.25 is Lx ~ 3 x 10*^ erg s~^ (within 
the measurement radius). We conclude that it is likely 
that this extended source is a cluster at z«0.25 with 
low Lx- 



9 DISCUSSION 

The evolution of cluster entropy was investigated by 
comparing the gas entropy measured at 0.1i?2oo in 
C1J 0046. 3+8530 w i th th at measured in local systems 
by IPonman et alJ l|2003r . Under self-similar models, 
entropy is simply proportional to temperature. How- 
ever with a large sample o f systems with high-quality 
data, IPonman et al.l l|2003|) found that this was not the 
case. To remove the effect of temperature dependence, 
C1J0046. 3+8530 was compared with a data bin of eight 
systems in the temperature range 3.6 — 4.6 keV with 
kT = 4.1 + 0.4 keV {z = 0.036). The entropy at 0.1ii2oo 
in C1J0046. 3+8530 is 5 = 180 ± 24 cm^ and the mean 
entropy of the eight local systems at the same radius is 
5 = 327 ±42 cm^ 

The mean density within a given overdensity ra- 



dius (relative to the critical density) is proportional to 
H{z)^, assuming self-similar scaling. Thus, the redshift 
dependence of scaled entropy is given by {HoE{z))~''''^, 
where 



E{z) 



(1 + 2:) 1 + zQm + 



Oa 



(1 + ^) 



fi/ 



1/2 



.(3) 



Using this formulation, the entropy in 
C1J0046. 3+8530 was scaled to the mean redshift 
of the local systems (assuming that the systems 
virialised at the redshift at which they are observed), 
giving S — 285 ± 40 cm'^, which is consistent with 
the local systems. This result is consistent with a 
simple self-similar evolution of entropy with redshift, 
a result also seen in the m assive, relaxed cluste r 
C1J1226.9+3332 at 2 = 0.89 JMauehan et al.ll2004^ . 
This result will be investigated further with a larger 
sample of high-z clusters in a future paper. 

In the context of other high-redshift {z > 0.6) clus- 
ters observed with Chandra and XMM - Newton (e.g. 
IVikhlinin et alJl2002l: lArnaud et al|l2002l: iHolden et al.l 



:200J: iMauehan et al.M2003l . l2004l: iRosati et alJl2004^ . 
C1J0046. 3+8530 has an average temperature and mass. 
The long exposure time of these two observations en- 
abled the system to be studied in more detail than 
is usually possible at high redshift. The measure- 
ment of the gas temperature structure and detection 
of emission to approaching the virial radius is rela- 
tively rare; temperature proflles have been measured 
in only three other z > 0.6 clusters, MS 1054-0321 
djoltcma et "8111200 1^ . RX J1120. 1+4318 ("Arnau d et ail 
2002) and C1J1226.9+3332 |||4auahan ct al. 200^J! 

The total gravitating mass of 3.0lo;5 x IO^'^Mq de- 
rived for the system assumed isothermality, which is 
consistent with the observed temperature proflle, and 
also supported by the hardness-ratio map. Thus this 
mass should be reasonably reliable. However, it is also 
possible that there is a temperature gradient present, 
which was not detected because of the effect of the 
telescope PSF. In a sample o f 66 local systems with 
measured temperature proflles ISanderson et al.l i2003l) 
found that the incorrect assumption of isothermality 
leads to an average overestimate of the total mass by 
« 30% at 7?2oo. Thus this is a reasonable estimate of 
the possible systematic error in our total mass measure- 
ment. 

The measured metal abundance in C1J0046. 3+8530 
is consistent with the value of O.3Z0 found in 
samples of low- and moderate- redshift clusters {e.g. 
iMushotzk v fc LoewensteinI |l993). In particular, the 
data impose a 3a lower limit oi Z > 0.17Zq, which sup- 
ports the growing evidence of a high redshift (a > 1) of 
enrichment of the ICM (e.Q. IMushotzkv fc LoewensteinI 
11991: iTozziet alJl2003^ . 
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The gas mass fraction of fgas = 0.16 ± 0.04 de- 
rived at -R200 for C1J0046. 3+8530 should be a reason- 
ably reliable measurement, because of the large extent 
of the surface-brightness and temperature profiles. This 
value is in good agreement with the mean gas mass 
fraction of fgas = 0.13 ± 0.01 (also at R200) found by 
ISanderson et alj i2003r in their large sample of local 
systems. In a sample of 6 relaxed clusters in the red- 
shift range 0.1 < 2 < 0.5 observed with Chandra , and 
with reliabl e mass es confirmed by gravitational lensing, 
I Allen et al.l (|2002r found a mean gas mass fraction of 
0.113 ± 0.005. While those measurements were made 
at a radius of r25oo (corresponding to a mean internal 
density of 2500 times the critical density), the value we 
derive for C 1J0046. 3-1-8530 at -R200 is consistent with the 
I Allen et al.l (|200^ measurements. 

In addition, the gas density slope parameter 
(/3) measured here for C1J0046. 3-1-8530 is consistent 
with those of s imilar temperature clusters in the 
ISanderson et alJ (|200o|) sample. Overall, these results 
suggest that the enrichment and distribution of the ICM 
is similar in at least some high-redshift clusters to that 
in local clusters. 



10 CONCLUSIONS 

A detailed XMM-Newton study of the high-redshift 
galaxy cluster C1J0046. 3-1-8530 (z=0.624) was con- 
ducted. The circular appearance of the isophotes, and 
the low ellipticity of the best-fitting surface bright- 
ness model suggest that the ICM is reasonably relaxed. 
A temperature profile and hardness-ratio map show 
no significant temperature structure. The metal abun- 
dance, gas-density profile slope, and gas- mass fraction of 
C1J0046. 3-1-8530 are all consistent with those measured 
in local clusters. The entropy of ClJ0046.3-f 8530 mea- 
sured at O.I-R200 is consistent with a self-similar model, 
in which high-z clusters are identical to their low-z coun- 
terparts, with their properties scaled by the increasing 
(with redshift) density of the universe in which they 
formed. 

The properties of C1J0046. 3-1-8530, along with 10 
other high-redshift clusters from the WARPS, will be 
used to investigate the evolution of the scaling relations 
between cluster properties in a forthcoming paper. 
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